Wide-field Fourier transform spectral imaging 
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We report experimental results of parallel measurement of spectral components of the light. The 
temporal fluctuations of an optical field mixed with a separate reference are recorded with a high 
throughput complementary metal oxide semi-conductor camera (1 Megapixel at 2 kHz framerate). 
A numerical Fourier transform of the time-domain recording enables wide-field coherent spectral 
imaging. Qualitative comparisons with frequency-domain wide-field laser Doppler imaging are pro- 
vided. 



Many coherent spectral detection schemes using a sin- 
gle detector (or balanced detection) to detect tempo- 
ral fluctuation spectra in an optical mixing configura- 
tion rely on Fourier Transform spectroscopy (FTS) for 
signal measurement [3, 0]. They provide a high spec- 
tral resolution and shot-noise sensitivity. They allow to 
shift away the 1// noise of laser intensity fluctuations 
since the measurement is done with GHz-bandwidth de- 
tectors. Most imaging configurations require a spatial 
scanning of the beam, but two approaches to parallel co- 
herent spectral imaging with a solid-state array detector 
were presented recently : full-field laser Doppler imag- 
ing (LDI) la, Jj, [5J and frequency-domain wide-field LDI 
(FDLDI) [a, LD, ls| • In the former approach, the tempo- 
ral fluctuations of an optical object field impinging on 
a complementary metal oxide semi-conductor (CMOS) 
camera are recorded. Spectral imaging is done by cal- 
culating the intensity-fluctuation spectrum by a Fourier 
transform (FT). One major weakness of this approach 
lies in its inapplicability in low-light conditions. The lat- 
ter approach uses a spatiotemporal heterodyne detection, 
which consists in recording an optical mix of the object 
field with an angularly tilted and frequency-shifted lo- 
cal oscillator (LO). It enables to measure spectral maps 
with a high sensitivity but requires to acquire the spectral 
components sequentially by sweeping the LO frequency. 
We present an alternative approach, designed to combine 
the advantages of both methods. It uses the properties of 
digital off-axis holography and FTS to enable exploring 
of the temporal frequency spectrum of the object field. 
Basically, the parallel spectral imaging instrument pre- 
sented here uses a CMOS camera to record the intensity 
fluctuations of an object field mixed with a separate ref- 
erence (LO); the field spectral components are calculated 
by FTS. 

The experimental setup is based on an optical interfer- 
ometer sketched in FigEQ A CW, 80 mW, A = 658 nm 
diode (Mitsubishi ML120G21) provides the main laser 
beam (field E^, angular frequency cJl)- A small part of 
this beam is split by a prism to form a reference (LO) 
beam, while the remaining part is expanded and illu- 
minates an object in reflection with an average incidence 
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FIG. 1: Setup. LD : single mode laser diode. M : mirror. BE 
: beam expander. BS : beam splitter. E : object field. Slo : 
local oscillator field. 



angle a « 45°. The object is made of a USAF 1951 target 
set in front of a 4 mm-thick transparent tank, filled with a 
non dilute intralipid (TM) 10% emulsion. To benefit from 
heterodyne gain, the field scattered by the object, E, is 
mixed with the LO field E LO ( \E LO \ 2 /\E\ 2 ~ 10 3 ), and 
is detected by a CMOS camera (LaVision HighSpeedStar 
4, 10 bit, 1024 x 1024 pixels at uj s /{^) = 2.0 kHz frame 
rate, pixel area d P i X 2 with d p i X = 17.5 /iin, set at a dis- 
tance d = 50 cm from the object. A 10 mm focal length 
lens is placed in the reference arm in order to create an 
off-axis (8 w 1° tilt angle) virtual point source in the 
object plane. This configuration constitutes a lensless 
Fourier holographic setup [9(. 

In the detector plane, the LO and object fields are: 



Eho(t) = £Loe 
E(x, y, t) = £(x, y, i)e ML * + c.c 



(1) 



where c.c. is the complex conjugate term. The LO beam 
is a spherical wave propagating along z, and thus the LO 
field envelope Elo does not depend on x, y, t. The object 
field envelope £, which contains information on the ob- 
ject shape, and which may exhibit speckle, depends on 
position x, y. It also depends on time t because of dy- 
namic scattering. The intensity I recorded by the camera 
can be expressed as a function of the complex fields : 



I(x,y,t) 



\E(x,y,t) + E^oiW 
= \£{x,y,t)\ 2 + \£ LO \ 2 



(2) 



2 



+E(x, y, t)£* LO + £*(x, y, t)£ LO 

where A is the time average of A over the optical pe- 
riod. The camera records the interference pattern of the 
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FIG. 2: Lensless off-axis holograms reconstructed in the tar- 
get plane. Images represent the field intensity, displayed in 
arbitrary logarithmic scale, (a) image obtained from a sin- 
gle hologram, (b) image obtained from the difference of two 
holograms. 

LO field £lo with the signal field £ , the recorded signal 
I(x, y, t) is the numerical hologram of the object that can 
be used to reconstruct the object image (loj . Because of 
the lensless Fourier holographic configuration, the recon- 
structed image field amplitude £ is obtained from / by a 
two-dimensional (2D) FT [H E2, : 



£{k x ,k y ) = FT 2£> I{x,y) 



(3) 



FigE] shows intensity images of the USAF target 
(i.e. | £(k x ,k y )\ 2 ) displayed in logarithmic scale (arbi- 
trary units). Fig[2k is obtained from a single frame I(ti) 
recorded at time t\. The USAF target is not visible 
because the noise is too large. To lower the noise, we 
have recorded two frames I{t\) and Ifa) at instants t\ 
and ti, and substracted them. By making the differ- 
ence of the two holograms, the noise components which 
do not vary with time (like the LO beam noise and the 
CMOS dark signal noise) cancel-out, whereas the Eq(2] 
holographic cross terms {££* LO and £*£lo) do not van- 
ish, because the signal field envelopes £{t\) and £(£2) are 
(at least partially) decorrelated in both amplitude and 
phase from one frame to another as a consequence of dy- 
namic backscattering by the intralipid emulsion. Fig[2j3 
shows the reconstructed intensity image (|£'| 2 ) obtained 
from the difference of two frames. One can notice that 
the last 4 terms of Eqf2] are visible on FigfSJa. The 
true image (white circle 1) corresponds to the cross term 
£££ , while the twin image (white circle 2) corresponds 
to £*£lo- Because of the lensless configuration, the true 
and twin images are on focus in the same reconstruction 
plane (i.e. the reciprocal plane of the detector). To pre- 
vent overlapping of the true, twin and zero order images 
in the off-axis holographic configuration, the true image 
size (circle 1 of diameter 0.75 cm 409 pixels) is ~ 2.5 x 
smaller than the total 1024 pixels field corresponding to 
Xd/dpi x = 1.87 cm [10]. This means that, on average, 
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FIG. 3: 128 x 256 pixels images of temporal frequency com- 
ponents of the object field, measured with the presented FTS 
instrument using a CMOS detector: (a) to (d), and the wide 
field FDLDI setup using a CCD detector: (f) to (i). (e) and 
(j): average over frequencies. Color axis is in logarithmic 
arbitrary units. 



one speckle grain is about 2.5 pixels. The light collection 
efficiency is 2.5 2 ~ 6x lower than with on-axis (or inline) 
holography or with homodyne detection (for which one 
speckle = 1 p ixel). \£lo\ 2 contributes to the zero order 
image [l^, EH • It yields the very bright region in the cen- 
ter (null spatial frequency) of FigfS^i and Figf2jD (arrow 

3) . Contrarily to \£lo\ 2 , the \£\ 2 term is not flat-field. 
It yields the broad spot in the center of Figf2b (circle 

4) . Because the brownian spectrum is narrower than the 
Nyquist frequency of the time-domain sampling, £ and 
/ vary not too fast in time to be sampled properly. It is 
then possible to record with the CMOS camera the time 
evolution of intensity fluctuations in time. From a se- 
quence of CMOS images, one can thus extract the Fourier 
temporal frequency components of the holographic sig- 
nal, and reconstruct images from these spectral compo- 
nents. We have recorded a data cube made of a sequence 
of TV = 2048 images at a framerate uis/(2ir) = 2 kHz. A 
3D numerical FT (2D for space, ID for time) was applied 
to this data to calculate spectral component maps of the 
object field envelope in the target plane : 



£(kx 1 ky 



FT 1D £(k x ,k y , t) = FT au I(x, y, t) (4) 



3D 



The FT along the temporal dimension is used to cal- 
culate spectral maps of the object field in quadrature 
(amplitude and phase), and the FT in the spatial di- 
mensions yields the field distribution in the object plane 
(image). Since we have performed a discrete FT, the 
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FIG. 4: Spectrum of the field dynamically backscattered by a 
non dilute suspension of intralipid 10% obtained by averaging 
the object field intensity over 50 x 50 pixels. FTS (solid 
gray curve) and FDLDI (points) spectra. Horizontal axis is 
the frequency u/(2n) in kHz, vertical axis is signal in linear 
arbitrary units. 



2048 frequency points u) are linearly spaced between the 

Nyquist frequencies ±1.0 kHz. The measurement time 

of the 1024 x 1024 x 2048 data cube is ~ 1 s, and the 
3D 

FT calculation time on a personal computer is about 
1 hour nowadays. Fig[3ja) to[3£d) show the images of the 

object field intensity \£\ 2 in the target plane for the fre- 
quency components to/(2ir) = (a), 187.5 (b), 500.0 (c) 
and 937.5 Hz (d) (128 x 256 pixels crops of the total holo- 
gram, displayed in logarithmic scale). For u> — (a), the 
LO beam noise is dominant and the target is not visible. 
For to 0, the USAF target is visible but the bright- 
ness and SNR of the image will decrease with frequency 
((b) to (d)). Fig[3]Je) shows the image obtained by aver- 
aging over all frequencies. We have compared these re- 
sults with wide-field FDLDI images 0, 0] obtained with 
a charge-coupled device (CCD) camera (PCO Pixelfly: 
1280 x 1024 pixels, framerate: 8 Hz) with four-phase de- 
modulation over 32 images per spectral point, in the same 
experiment. FigOshows the 128 x 256 pixels FDLDI im- 
ages at (f), 187.5 (g), 500.0 (h) and 937.5 Hz (i), while 
image (j) corresponds to the average over all frequen- 
cies. The USAF target is seen on all the images. For 
lu = 0, the target appears as a contrast-reversed image 
6f . Since the pixel size of the CCD camera (6.7x6.7 ^im) 
is smaller than its CMOS counterpart (17.5/zmx 17.5/im), 
the extension of FDLDI image is larger (the acceptance 
angle of the receiver is proportional to the inverse of the 
pixel size). The white dashed rectangle of Fig[3j corre- 
sponds to the CMOS-imager field of view. Although the 
number of recorded spectral points was kept low for the 
FDLDI measurement compared to the FTS scheme (64 
vs. 2048), the total measurement time was much greater 
(256 seconds vs. 1 second). This difference is due to the 
throughput discrepancy between the CCD (1.3 Mpixel @ 
8 Hz) and the CMOS (1.0 Mpixel @ 2 kHz) receivers. 

We have computed the frequency spectrum of the light 
diffused by the intralipid emulsion with FTS. This spec- 
trum is obtained by averaging the object field intensity 



over a 50 x 50 pixels region of the reconstructed image. 
The lineshape is plotted on Fig[4] as a solid gray curve. 
We have compared its shape with the one obtained with 
the FDLDI technique (Fig|4] points) . The agreement is 
good except in the tails of the spectrum. The FTS fre- 
quency response is imperfectly fiat, because of the CCD 
finite ex pos ure time (1/tus) that yields signal low pass 
filtering lH, 1(3]. Moreover, because the signal temporal 
evolution is sampled at 2 kHz, temporal sampling aliases 
and spectrum overlap are expected around the Nyquist 
frequencies ±1 kHz. 

In this Letter, we have shown that the spatiotemporal 
heterodyne detection recently introduced @, S @| can be 
adapted to a wide-field Fourier transform spectral imag- 
ing scheme with a high throughput array detector. By 
using an off-axis optical mixing configuration, the object- 
LO fields cross terms are shifted away from center of the 
detector reciprocal plane (k-space) , contrarily to the ob- 
ject and LO self-beating contributions, which remain un- 
shifted. It is then possible to reject the local oscillator 
and the object field self-beating contributions accounting 
for noise. The heterodyne gain provided by optical am- 
plification of the object field by the LO field is essential 
for a high frame rate camera measurement in low-light 
conditions, since the object field intensity decreases with 
the camera exposure time. The ability to filter-off the 
LO beam noise, yields an optimal sensitivity of 1 photo- 
electron of noise per pixel. This limit has been reached 
with 4-phase detection (l7J, which consists of a discrete 
Fourier transform on 4 data points to calculate a sin- 
gle frequency component of the object field. Here, the 
expected noise limit is the same for each frequency com- 
ponent of the object field obtained by discrete Fourier 
transform. This method might find applications in dy- 
namic light scattering analysis of colloidal suspensions 
and microfluidic systems. 
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